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Nomenclature 


A cross-sectional area of regenerator, 

Aign cross-sectional area of regenerator, 
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molecular weight of Ha, kg kmoT^ 

pressure in regenerator, atm 

pressure in riser, atm 

rate of the i* reaction, i=l to 9 (riser); i = 10 to 12 (regenerator), 
kmol (kg catalyst)’’ s’’ or kmol m’^s’’ 
universal gas constant, J K’’ kmol’’ 


temperature of air fed to the regenerator, K 

base temperature for heat balance calculations, K (assumed, 866.6 K) 


TriSjtop 


Greek 


temperature of dilute phase at any location, K 
temperature of gas oil feed, K 
temperature (uniform) of dense bed, K 
temperature of riser at any location, K 
temperature at top of riser, K 
temperature of spent catalyst [= Tris,top - ATst], K 
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(kg catalyst entrained in dilute phase) (kg fluidizing vapor)'’ 
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height of the dense bed, m 
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total height of the regenerator, m 

stoichiometric coeflBcient of j* species in i”^ reaction, based on mass 
CO/CO 2 ratio at catalyst surface in regenerator (Table 2) 
void fiaction in riser or regenerator at any location 
void fiaction in the dilute phase at any location 
density of solid catalyst (not including void fiaction), kg m'^ 



Pdil 

*3 

density of catalyst in the dilute phase, kg m' 

Pf 

density of fluidization vapor, lb 

Pg 

density of gas phase in the regenerator, kmol 

Pp 

density of catalyst particle (solid), lb ft ^ 

pv 

density of vapor at any location, kg m 

<!> 

activity of the catalyst 

Subscripts 

ij 

i*^ or lump; 1 ; gas oil, 2; gasoline, 3: LPG, 4; dry gas, 5; coke 


ABSTRACT 


This study provides insights into the optimal operation of one of the most important 
refinery units, namely, the fluidized-bed catalytic cracking unit (FCCU). A five-lump 
model is used to characterize the feed and the products. The model is first timed using 
some industrial data. An elitist non-dominated sorting genetic algorithm (NSGA-II) is 
used to solve a few meaningful multi-objective function optimization problems. It is 
observed that the non-dominating Pareto set of optimal solutions can be broadly 
classified into two regions: the partial-combustion and the full-combustion problems. A 
three-objective function optimization problem is solved. The objective functions used 
are: maximization of the gasoline yield, minimization of the air flow rate and 
minimization of the percent CO in the flue gas. The flow rate of the gas oil feed to the 
industrial FCCU is taken as 29 kg/s. The decision variables and several important state 
variables corresponding to the optimal conditions of operation, are obtained. The 
procedure used is quite general and can be applied to other industrial FCC units. The 
optimal results obtained here provide physical insights that can help one in obtaining 
and interpreting such solutions. 



CHAPTER 1 


INTRODUCTION 


Increase in the demand for gasoline, LPG and diesel over the last several 
decades have led to major improvements in refinery operations. The fluidized-bed 
catalytic cracking unit (FCCU) is an important conversion unit in most integrated 
refineries. Several studies have been reported in the open literature that deal with 
various aspects of FCCU. These include their modeling, simulation, kinetics, on-line 
optimization and control. Avidan and Shinnar* reviewed the developments and 
commercialization of catalytic cracking in detail. Different workers have discussed the 
kinetics^"^ in the reactor and the regenerator and have modeled’"^ these units separately, 
while a few^°'^^ have developed an integrated model for the reactor-regenerator system. 

In case of the reactor, various models are available that describe the feed and the 
products in terms of different kinetic Tumps’. These include describing the products in 
terms of two^^, three^, four'* or five^’^ lumps, while considering the feed as a single 
lump. Researchers'^'^' have obtained some useful results using these models. The need 
to evaluate the dry gas, LPG, and the gasoline yields led to an increase in the number of 
lumps from two to five. Jacob et al."* proposed a ten-lump kinetic scheme that 
characterizes the feed in terms of eight different lumps. The product lumps include, in 
addition, gasoline and (coke + gas) lumps. Cerqueira et al.^^ proposed a twelve-lump 
kinetic scheme with coke and gas as separate lumps. The advantage of ten- and twelve- 
lump kinetic schemes is that the rate constants are independent of the feed composition, 
but the main problem is that a relatively large number of kinetic parameters need to be 
evaluated using e>q)erimental or industrial data that are not easily available. 



In the FCCU, the catalyst loses its activity due to the deposition of coke. 
Voorhies Nace et al.^'* and John et al.^^ related the catalyst activity to the time-on- 
stream, Froment and Bischofi?^ have related this decay to the amount of coke deposited 
on the catalyst, while Beeckman and Froment^^ have accounted for the pore size 
distribution and the arrangement of the pores in the catalyst. 

Several models are available in the literature for the regenerator, e.g., simple 
models,’’® grid-effect models,* two-region models’*’^^ and the bubbling bed model.^° 
Several studies^^'^^’^^-^^ have been reported on the integrated models that include both 
the reactor and the regenerator. It is clear that a considerable amount of effort has been 
put into the modeling of FCCU. A large number of complex equations need to be 
solved simultaneously to obtain results numerically. 

Several studies have been carried out on the optimization of FCCUs. Dynamic 
optimization of an FCCU was carried out by Davis et al.^' and Webb et al.'" Rhemann 
et al. carried out online control and optimization of FCCUs. McFarlane and Bacon 
used the adaptive optimizing control strategy to study the dynamic history of 
input/output variables with respect to changes in the ambient temperature and the ‘coke 
factor’. They used a profit function as the objective function in their study and the 
regenerator temperature, reactor temperature, catalyst circulation rate and the air supply 
rate as the four decision (optunizing) variables. Chitnis and Corripio^’ used SMCO 
(supervisory multivariable constrained optimization) to optimize FCCUs. The aim of 
their study was to minimize the objective function that included a cost term, a 
constraint term and a ‘move-suppression’ term. They foEowed the study done by 
McFarlane and Bacon. Ellis et al.^‘ used SQP (successive quadratic programming) to 



study the sensitivity of the economic objective fiinction for online optimization. They 
used the regenerator temperature, riser temperature, gas oil feed flow rate and stack gas 
oxygen concentration as the decision variables. Khandalekar and Riggs^^ used the 
optimization method of Nelder and Mead^^ for the inner cycle to get the regenerator and 
riser temperatures, and the Fibonnaci method"^*^ for the outer cycle to obtain the optimal 
feed rate. The decision variables used by Khandalekar and Riggs were the same as 
those used by Ellis et al. Recently, Ramasubramanian et al."*^ used the optimization 
procedure of Luus and Jaakola"^^ to maximize a performance index defined by them. 
They used nine different decision variables in their study. In addition to the decision 
variables used in the previous studies,^^’^^’^^ they used the wet gas compressor valve 
position, reactor-regenerator pressure difference, governor set-point for lift air blower, 
combustion air blower suction valve position, and the spill air valve position as 
additional decision variables. More recently, Zhao et al.'*^ used artificial neural network 
(ANN) and modified genetic algorithm (a modification of the conventional to 

optimize the operation of FCCUs. They used ANN for describing the relation between 
conversion ratio and the influencing factors. They found that ANN is better than 
available empirical models. Till now, all the optimization studies of FCCUs have used 
only a single objective function. In thus study, however, we use an adaptation of 
genetic algorithm (NSGA-II), as developed by Deb and coworkers"^®’^^ to optimize 
industrial FCCUs using more than one objective function. Such studies, we believe, are 
more realistic. Indeed, multi-objective optimization has been carried out for several 
industrial units"^^'^^ in chemical engineering using an earlier adaptation of GA 
developed by Srinivas and Deb'^^ (NSGA-I). 



This is the &st study in the field of chemical engineering in which we are 
using NSGA-II,'^®’^^ a more recent and improved version of NSGA-I,"^^ to obtain 
solutions for the FCCU, which is well known to be a computationally intensive system. 
Also, this is the first attempt in the field of FCCU system optimization where more than 
a single objective function has been used so as to obtain more meaningful and realistic 
results. 

Process Description 

The FCCU consists of two major units, a reactor/riser and a regenerator, as 
shown in Figure 1. The feed to this unit, consisting mainly of heavy gas oil, vacuum 
gas oil or coker oil, is preheated to a temperature, Tjbcd, and gets mixed with the hot 
regenerated catalyst at temperature, Ti^, at the base of the riser. The feed is mixed with 
a small amount of steam, which helps in good atomization and reduction in coke 
formation by reducing the partial pressure of hydrocarbon vapors. All cracking 
reactions taking place in the riser are endothermic. Lighter hydrocarbons are produced 
in the riser. Byproduct coke is formed simultaneously which deposits on the catalyst 
surface and reduces the activity of the catalyst. The amount of coke formed depends 
upon the coking characteristics of the feedstock and also on the operating conditions. 
At the top of the riser, the product vapors are separated from the catalyst in the 
disengaging section (separator in Fig. 1) of the reactor and sent to the main fractionator 
for separation. The hydrocarbon left on the catalyst surface is removed using steam 
stripping. The stripped spent catalyst at temperature, Tsc, is sent to the regenerator. The 
catalyst is fluidized here using air at temperature, Tair. The hydrogen and carbon on the 
regenerator catalyst react with oxygen in the air, resulting in the formation of carbon 



Figure 1 Schematic diagram of fluid catalytic cracking unit 


To main 




monoxide, carbon dioxide and water. The combustion reaction is exothermic. The 
height of the fluidized bed (dense bed) is Zbed- There is a catalyst disengaging space 
(Zda) above the dense fluidized bed. This is referred to as the dilute bed. Entrained 
catalyst is finally separated from the gas in the cyclones at the top of the regenerator. 
The exothermic heat associated with the burning of the coke is used to keep the catalyst 
at high temperatures. This supplies the heat required for the vaporization and cracking 
of the feed in the riser. Losses in catalyst activity due to coking are temporary. 
However, exposure of the catalyst to high temperatures and steam can result in its 
permanent degradation. Attrition results in the loss of catalyst with the flue gas, and is 
made up for using a supply of fresh catalyst. 

Model of the FCCU 

In the present study, the reactor model of Dave^^ is used because of its 
simplicity and applicability to industrial units. The five-lump kinetic scheme of 
Ancheyta et al.^ is used. This kinetic scheme is simple but somewhat inadequate since 
several of the rate constants depend on the feed composition. In contrast, the ten-lump 
model of Jacob et al."* does not suffer from this drawback, but is not used here since 
sufficient industrial data are not available in the open literature to ‘tune’ the several rate 
constants associated with it. Dave modified the original model of Ancheyta et al.^ by 
assuming that gasoline and LPG also convert to coke. This modified kinetic scheme is 
shown in Figure 2. The rate constants and the heats of reaction are given in Table 

I 14,32,57,58 

In the study of Dave, the FCC reactor/riser has been modeled as a steady state 
plug flow reactor with the assumption that gas oil cracking follows second order 




TABLE 1 


Kinetic and Thermodynamic Parameters'^ used for Reactor Modeling 


Rate constants 

Frequency 

factor*^^’^* 

Activation 

energy^’ 

(kJ/kmol) 

Heat of 
reaction^'* 
(kJ/kmol) 

ki 

14054.59^ 

57540 

45000 

k2 

2293.00^ 

52500 

159315 

ks 

390.95^ 

49560 

159315 

k4 

29.90^ 

31920 

159315 

ks 

65.40** 

73500 

42420 

k6 

0.00** 

45360 

42420 

k? 

0.00** 

66780 

42420 

ks 

0.32** 

39900 

2100^ 

k9 

0.19** 

31500'^ 

2100^ 


* (kg catalyst)'’ (kmol gas oil)'’s'' for ki - k 4 

m^(kg catalyst)'’s'’ for ks - k 9 
+ from Ref 32 


** from Refs. 32 and 58 




kinetics, while gasoline and LPG cracking reactions follow first order kinetics. The 
catalytic deactivation function for the cracking of vacuum gas oil is that provided by 
Yingxun.^^ Dave assumed the temperature drop aroimd the stripper to be about 10 °C, 
and adjusted the hydrogen-to-carbon ratio in the catalyst at the entry point of the 
regenerator, in order to match industrial data on the regenerator temperature. The 
regenerator is modeled using the two-region model proposed by Krishna and Parkin^ 
with some modifications.^^ Gas is assumed to be in plug flow throughout the 
regenerator bed and in thermal equilibrium with the surrounding solids. The catalyst is 
assumed to be well mixed in the dense bed and in plug flow in the dilute phase. 

The main reactions taking place in the regenerator are; 


c+\o,- 

Jh^CO 

(a) 

2 


C + O 2 


(b) 

C0 + -0, 
2 ^ 

-^co^ 

(heterogeneous CO combustion) (c) 

C0 + -0,- 
2 " 

-^co. 

(homogeneous CO combustion)(d) (1) 


The balance equations for O 2 , CO, CO 2 , as well as for H 2 O (reaction not listed 
in Eq. 1) are given in Appendix 1. The kinetic parameters used for the regenerator 
model are given in Table 2.^‘^ The complete set of equations^^ are summarized in 
Appendix 1. Table 3 gives the thermodynamic and other parameters while Table 4 
gives the design details of the industrial FCCU used in this work. These are slightly 
modified from those corresponding to an existing industrial unit,^^ for proprietary 


reasons. 



TABLE 2 

Kinetic Parameters used for Regenerator Modeling 


Parameter 

Frequency factor 

Activation energy 
(E/R, K) 

Pc* 

2512 

6795 

ko atm‘’s'’ 

1.0694x10** 

18890 

kl2c 

(kmol CO)(kg caty^s'^ m'^ 

117 

13890 

kl2h 

(kmol CO)m'^atm'^s'^ 

S.OTxlO^'* 

35555 


♦:P..-^=|3ooexp[-^] 
Ki RT 

— E 

kc = kio+ kn = kco exp [-^] 

RT 







TABLE 3 

Thermodynamic and Other Parameters'^ used for Simulation 


Parameter 

Value 

Cpo kJ kg'^ K'^ 

1.003 

Cp^kJkg-'K-' 

3.430 

Q^^^kJkg’ K’’ 

3.390 

kJkg-'K-' 

30.530 

kJkg-'K-' 



32.280 

Cp^.okJkg-^K-^ 

36.932 

Cpco kJ kg'^ K'^ 

30.850 

Cpco^ kJkg-'K-' 

47.400 

AHe^kJkg-^ 

350.0 

Hco kJ kmof^ 

1.078x 10^ 

Hco^ kJ kmoT* 

3.933 X 10^ 

kJ kmof' 

2.42 X 10^ 




TABLE 3 (contd. . .a) 


Parameter 

Value 

pokg m'^ 

1089.0 

CH(kgH2)(kg coke)"’ 

0.15 

MWgtisoil 

350 

M gasoline 

170 

MWi^PQ 

65 

IVlW (Jiy £t^g 

30 

MWcoke 

12 

♦ 

Otgas oil, gasoline 

2.06 

Ofgasoa,LPG 

5.39 

Otgas oil, diy gas 

11.67 

oil, coke 

29.17 

OCgasoline, LPG 

2.62 

OCgasoline^ diy gas 

5.67 








TABLE3(contd...b) 


Parameter 

Value 

Otgasoline, cjoke 

14.17 

J^G, dry gjts 

2.17 

OtLPG, coke 

5.42 

Dp ft 

2.0 X lO"* 

Xpt 

0.10 


* : Oij = forthwith liunp corresponding to the reaction from lump i to lump j (in 
mass tertns) 




TABLE 4 


Design Data for the FCCU Used^^ in This Study 


Parameter 

Value 

Riser length (m) 

37.0 

Riser diameter (m) 

0.685 

Regenerator length (m) 

19.4 

Regenerator diameter (m) 

4.5 

Inventory of catalyst in regenerator (kg) 

34,000.0 

Feed rate (kg/s) 

29.0 

Riser pressure (kPa) 

253.85 

Regenerator pressure (kPa) 

267.23 




The regenerator model of the FCCU is coupled with the reactor model, and for 
any specified operating conditions, it can predict the product yields, coke on 
regenerated catalyst, temperature of the regenerated catalyst, and the flue gas 
composition. The Runge-Kutta-Gill method^® is used to solve the differential equations. 
Guess values are assumed for the coke on regenerated catalyst, C,^c, and the 
temperature, Trgn, of the regenerated catalyst, and the equations are solved to give new 
values (iterates) for these two variables. Convergence is obtained using the Newton- 
Raphson method.^ A tolerance of 1.5°C is used for the convergence of Trgn, while a 
tolerance of 1 x 10"^ kg coke/kg catalyst is used in the case of Crgc. It was found that 
sometimes the model does not converge, particularly when the decision/input variables 
generated by GA (that uses a random number generator for this purpose) are practically 
infeasible. 

The model is tuned in order to describe the operation of the industrial FCCU 
more realistically. It was observed that the temperature at the base of the riser in the 
equation was always higher than the industrial values. Similarly, the model-predicted 
values of the dense bed temperature, T,gn, differed slightly from industrial values. In 
order to minimize these two deviations, empirical heat loss terms were introduced for 
the riser base (dotted box shown in Fig. 1) as well as for the dense bed of the 
regenerator. The appropriate equations are also included in Appendix 1. Introduction of 
these losses led to good agreement of model-results with four sets of industrial data.^^ 
Results for all four sets of data are shown in Appendix II. A similar loss term was 
incorporated by McFarlane et al.*^ and Han et al.^^, but only for the dense bed heat 


balance. 



Optimization 

In most real life problems, one needs to optimize several (incompatible) 
objectives functions. This could result in a set of several non-dominated (equally good) 
Pareto-optimaf ^ solutions, rather than a unique optimal point. The advantage of using 
multi-objective evolutionary techniques is that one can obtain the entire Pareto set in a 
single application of the algorithm. GA, as first developed, has been adapted to 
solve multi-objective problems. One such adaptation for problems involving decision 
variables that are single values, is the non-dominated sorting genetic algorithm (NSGA- 
1) developed by Srinivas and Deb.'*^’^^ Mitra et al.'**^ extended this technique to apply to 
problems where the decision (optimizing) variables are Junctions of time (trajectory- 
optimization problems). In the last few years, several studies have been reported on the 
multi-objective optimization of industrial units using NSGA-I. Examples in chemical 
engineering include the optimization of polymerization reactors,"*^’^*^ steam reformers, 
cyclone separators,^^ membrane separation modules,^^ etc. These have been reviewed 
recently by Bhaskar et al.^^^ Unfortunately, NSGA-I has been observed to foil in a few 
cases, e.g., PET reactors,^®’^® described by a complex set of equations. A 
modification'*®’^^ of this algorithm, the elitist non-dominated sorting genetic algorithm 
(referred to, here, as NSGA-II) is now available in the literature that incorporates the 
concept of elitism in order to make the algorithm more powerful. Details of NSGA-II 
are provided in Appendix III (it is assumed that the reader is aware of the preliminary 
details'*^ on GA and NSGA-I). 



CHAPTER! 


FORMULATION 


Several objective functions can be considered in any optimization study of 
FCCU in order to maximize its profitability and satisfy real-life constraints. In this unit, 
an increase in the yields of gasoline, LPG or diesel always leads to an increase in the 
profit. So, it is clear that the maximization of these yields should be taken as the 
objective function. Interestingly, an increase in the gasoline yield usually leads to an 
increase in the yield of LPG, and of dry gas, and so the use of gasoline yield as one of 
the objective functions is sufficient. The kinetic model^^ used here is not detailed 
enough to predict the yield of diesel, and so the latter is not included. In contrast, an 
increase in the yield of gasoline has an adverse effect on (increase in) the coke 
formation. This coke decreases the activity of the catalyst, and needs to be burnt off in 
the regenerator, requiring higher amounts of feed air. Burning of the coke results in the 
formation of CO and CO2. The former needs to be converted to CO2 in the dilute phase 
before the gases enter the cyclones. This after-burning of CO can produce very high 
temperatures. There are usually two options available: one is to carry out full 
combustion in the regenerator, so that the CO emitted in the flue gas is very low. The 
other is to carry out only partial combustion in the regenerator and allow emission of 
higher amounts of CO firom this unit. The flue gas is then sent to a CO reboiler where it 
is converted to CO2 before being emitted finally to the atmosphere. The full combustion 
mode of the FCCU requires large amounts of air supply. This automatically increases 
the operating costs. In addition, full combustion leads to much higher amounts of 




exothermic heat produced, which could create problems. There is also a need to bum 
off the maximum amount of the coke in the regenerator so as to keep the catalyst 
activity in the riser at high levels. Based on the above discussion, we could select three 
objective functions for the study. The first is to maximize the gasoline yield 
(profitability), the second is to minimize the CO in the flue gas (pollution limit) and the 
third is to minimize the air flow rate (operating costs). Consideration of all three 
objective functions simultaneously at the start of any multi-objective study is difficult 
to attempt and understand, and so two simplified, two-objective function problems are 
first studied. 

The mathematical representations of the two-objective function optimization 


problems for the FCCU are: 

Problem Nos. 1 and 2 : 

Max Ii (u) = gasoline yield (a) 

Min I2(U) = Fair (b) 

subject to (s. t.): 

700K<Tiga <950K (c) 

C^c<l% (d) 

CO in flue gas < 8 % (Problem 1), or 

CO in flue gas < 1000 ppm (Problem 2) (e) 

where 

U — {Tfeed, Tair, Fcatj Fair} (0 (2) 

The following bounds are used for the four decision variables 

575^Tfeed<670K (a) 



450<Tair <525K 

(b) 

115 < Feat <290 kg/s 

(c) 

11 < Fair <46kg/s 

(d) (3) 


The feed preheat temperature, Tfeed, is selected as a decision variable because it 
plays a major role in controlling the heat balance in the FCCU. The lower limit of Tfeed 
is selected so as to supply sufficient heat required for the cracking and vaporization of 
the feed. Too low a value of Tfeed results in a decrease in the riser temperature which, 
ultimately, decreases the yields of gasoline, LPG, etc. The upper limit’'^ of Tfeed is 
decided so as to prevent coking of the heating coils in the preheater. The riser top and 
the regenerator temperatures are functions of the catalyst circulation rate. Feat, and the 
air flow rate. Fair- The lower and upper bounds of the catalyst flow rate, Feat, are taken 
such that the catalyst-to-feed flow ratio, Fcat/Ffeed, (keeping Fau/Ffeed constant) lies 
between^'^ 4 and 10, respectively. Similarly, the bounds on the air flow rate, Fat, are 
taken such that the air-to-feed flow ratio, Faii/Ffeed, (keeping Fcat/Ffeed constant) lies 
between^"^ 0.4 to 1.6. Below and above these ratios, trivial steady states^"^ exist, that 
have no relevance in industrial operations. The lower and upper bounds of the air 
preheat temperature. Tan, are determined by the amount of heat that air can pick up in 
actual operations. Too low a value of Tan- results in a cooling of the regenerator, while 
too high a value has the opposite effect. The first end-point constraint (Eq. 2c) on Tiga is 
incorporated so as to maintain reasonable temperatures in the dense bed. Too low a 
value of T;tgn decreases the temperature of the regenerated catalyst, which causes a 
decrease in the riser temperature, which, m turn, leads to a decrease m the yields of 
gasoline, LPG, etc. Too high a value of T,ga reduces the catalyst activity rapidly and 



also causes damage to the catalyst in the regenerator. The second end-point constraint 
(Eq. 2d) is on the coke on regenerated catalyst. Too high a coke content on the 
regenerated catalyst results in a decrease in the catalyst activity, which automatically 
decreases the yields of gasoline, LPG, etc. The third end-point constraint (both 
possibilities in Eq. 2e) is put on the amount of CO in the flue gas. The two-objective 
optimization problem is solved in order to get the optimal conditions for operating the 
FCCU in the partial combustion mode (Problem 1) or the full combustion mode 
(Problem 2). The value used for CO in Problem 1 is somewhat higher than the 
industrial value available to us. 

It is found (see Results and Discussion later) that the maximum amount of CO 
allowed in the flue gas (Eq. 2e) plays an important role in determining the optimal 
conditions for FCCU operation. Too low a limit results ia large amounts of heat 
liberated m the regenerator (due to after-burning). This ultimately results in very large 
temperature increases in the dilute phase of the regenerator. This could cause damage to 
the cyclones. 

In order to develop further insight, we studied yet another two-objective 
function optimization problem: 


Problem No. 3 


Max Ii (u) = gasoline yield 

(a) 

Min I 2 (u) = % CO in the flue gas 

(b) 

700 K<T^ < 950K 

(c) 

Crgc<l% 

(d) 

where 




« = {Tfeed, Tair, Feat, Fair} 


(e) (4) 


The bovmds on the decision variables (Eqs. 3 a-d) are the same as in Problems 1 and 2 . 
The results of this study show that the solutions in the partial-combustion zone (with 
higher CO levels) are not picked up by NSGA-II since they are inferior to those in the 
complete-combustion zone (with CO at ppm levels). This is because one of the 
objective functions is to minimize the CO in the flue gas. 

In order to encompass Problems 1-3 into a single interesting multi-objective 
optimization problem, we solved the following three-objective function optimization 
problem: 


Problem No. 4: 

Max Ii (u) = gasoline yield 

(a) 

Min h (u) = Fair 

(b) 

Min I 3 (u) = % CO in the flue gas 

(c) 

s.t. 

700K<Tr^^950K 

(d) 

Crgc<l% 

(e) (5) 


The bounds on the decision variables (Eqs 3 a-d) are the same as in Problems 1-3. 

Most workers in the field have used the penalty-function approach while solving 
optimization problems with end-point constraints.^'^ This works when the process is not 
too complex. In the case of the FCCU, the model is quite complex and a different 
approach is being used. Since GA generates solutions randomly, the moment a solution 
is generated, it is checked to see if it is able to satisfy all the end-point constraints. If it 



does not, we do not accept it as a member of the population, and generate another 
chromosome till these conditions are satisfied. This enables all the members in a 
population to satisfy all the constraints. A similar thing is done when daughter 
chromosomes are obtained after crossover and mutation. This technique was found to 
be superior to the conventional penalty function approach'*® (penalty proportional to the 
deviation) for the present problem. It enables the generation of a wider variety 
(diversity) of solutions. However, it is a little slow in the beginning. 



CHAPTERS 


RESULTS AND DISCUSSION 

A computer code was written in MS Fortran Power Station^'^ to obtain the results 
of the multi-objective optimization problems. The code was tested (using a single 
chromosome) for simulation of a few runs for which industrial data was available to us, 
and compared with results obtained by Dave,^^ The parameters describing the FCCU are 
given in Table 4. The initial guesses of Tiga and C,gc were 900 K and 0.002 kg coke/kg 
catalyst, respectively. The parameters associated with GA for Problems 1-4 are given in 
Table 5. The computational time required for solving Problem 1 is 12 hr on a Pentium 4 
(1.7 GHz) machine. The results for Problems 1 and 2 are shown in Fig. 3. If we go from 
any one point to another in Fig. 3 a, we find that one objective function, e.g., gasoline 
yield, improves (increases) while the other, the air flow rate, worsens (increases). These 
plots are, therefore, Pareto sets. The decision variables associated with each of these non- 
dominated solutions are shown in Figs. 3 b, c, and d. Figs. 3 e-q show several important 
state variables under these conditions. 

We first discuss the results iox Problem 1. It can be observed from Figs. 3 b-d 
that the air and the feed preheat temperatures for most of the optimal solutions lie at their 
upper bounds, and that the catalyst flow rate is the most sensitive of the decision 
variables. A close study of Fig. 3 d shows that the catalyst flow rate first increases to its 
upper boimd (with simultaneous changes in the air and feed preheat temperatures) as the 
air flow rate increases, and then drops sharply, before increasing again. In this second 
region, the other two decision variables remain essentially constant at their upper bounds. 
Some scatter is observed in Fig. 3 b. This is because of the relative insensitivity of the air 



TABLES 

GA Parameters'^^ Used in this Study 


Parameter 

Problems 1, 

2, and 4 

Problem 3 

Total chromosome length (Ndr) 

40 

40 

Number of generations (Nga) 

50 

90 

Population size (Np) 

50 

100 

Crossover probability (po) 

0.95 

0.95 

Mutation probability (pm) 

0.05 

0.05 
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Figure 3 Pareto set, decision variables and state variables of optimal solutions 
for Problem 1 (unfilled circles) and Problem 2 (filled circles) 
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Figure 3 (contd.) Pareto set, decison variables and state variables of optimal solutions 
for Problem 1 (unfilled circles) and Problem 2 (filled circles) 
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Figure 3 (contd.) Pareto set, decision variables and state variables of optimal solutions 
for Problem 1 (unfilled circles) and Problem 2 (filled circles) 






preheat temperature in influencing the Pareto. Similar insensitivities of decision variables 
leading to scatter have been observed in our earlier studies, too, and one could, if one 
wished to, easily select constant values of such variables and solve the multi-objective 
optimization problem again. Increase in the feed preheat temperature results in an 
expected increase iti the temperature at the base of the riser (Fig. 3 f). An increase in the 
air preheat temperature, however, increases the heat input to the regenerator which causes 
a rise in the temperature of the regenerated catalyst (Fig. 3 n) and, ultimately, increases 
the riser base temperature. Similarly, an increase in the catalyst flow rate results in an 
increase in the cracking reactions. This leads to an increase in the total conversion (Fig. 3 
e) as well as in the individual yields of gasoline, dry gas and LPG (Figs. 3 a, j, k). 
Simultaneously, this leads to an increase in the amount of coke formed (Fig. 3 i) in the 
riser. In order to satisfy the constraint on Ci^c (Eq. 2 d), the decision variables at lower 
values of the air flow rate in the Pareto set (Fig. 3 a) take on lower values compared to 
values at higher air flow rates. Increase in any of the decision variables (except the ah 
flow rate) in this region (lower ah flow rates) of the Pareto results in higher coke 
formation in the riser. Under these conditions, the amount of oxygen available both at the 
feed-end and at the top of the dense bed (Fig. 3 p) in the regenerator, is low. This results 
in the coke on the regenerated catalyst (Fig. 3 1) to lie at its maximum permissible value. 
Also, this reduces the after-burning of CO in the dilute phase of the regenerator, and 
leads to lower values of the temperature increase (Fig. 3 g) there. In contrast, at higher ah 
flow rates on the Pareto set, more oxygen becomes available at the end of the dense bed 
(Fig. 3 p). This causes after-burning of CO and raises the temperature in the dilute phase. 
This also leads to a decrease in the coke on the regenerated catalyst, even though more 
coke is formed in the riser under these conditions. It is observed from Fig. 3 h that the CO 



in the flue gas is at reasonably high levels of around 4 - 8% for all the optimal solutions 
in Problem 1 (partial-combustion problem). 

Problem 2 constraints the CO in the flue gas to lie at low levels (complete- 
combustion problem). Fig. 3 a shows the Pareto set for this problem too. Clearly, higher 
air flow rates are necessary to achieve this. The corresponding decision and state 
variables for these solutions are shown in Figs. 3 b-q. All the parameters (conversion, 
gasoline yield, coke formed in the riser, dry gas yield, LPG yield, temperature at the base 
of the riser) show almost similar trends as for Problem 1. However, there is a dramatic 
increase in the temperature in the dilute phase (Fig. 3 g). This is because the air supplied 
is quite high compared to values in Problem 1, resulting in the availability of more 
oxygen at the end of the dense bed (Fig. 3 p) in the regenerator for the after-burning of 
CO in the dilute phase. The amount of CO in the flue gas now reduces to ppm levels 
(full-combustion). Also, coke on the regenerated catalyst is at negligible levels under 
these conditions, since almost all the coke is burnt. 

It may be argued that the non-dominating solutions of Problem 2 should form a 
subset of Problem 1. This is not so. Fig. 3a shows that for the same gasoHne yield. 
Problem 1 leads to lower flow rates. Since the air flow rate is minimized, the algorithm 
does not pick up the points of Problem 2 while solving Problem 1 (unless it is forced to 
by the constraint on CO in the flue gas). The computational time required in Problem 2 is 
extremely high (for 10 generations, it takes over 96 hrs) because the algorithm has to re- 
generate new chromosomes whenever all the end-point constraints are not satisfied, and 
in this problem, this happens very frequently. 



As mentioned earlier, the solutions in the partial-combustion zone are not picked 
up by the optimization algorithm in Problem 3 (refer Appendix IV), since they are 
inferior to those in the complete-combustion zone (because the CO is being minimized). 

With the physical insights (and computational skills to obtain the correct Pareto 
solutions) now developed, we proceed to solve the more difficult Problem 4 (Eq. 5) 
involving three-objective functions. The population size was increased to 100 for this 
problem, and the number of generations was taken as 90. The computer time was about 
48 hr. The code for the three-objective function optimization problem was re-run with 
two of the objective functions identical (so as to reduce to Problem 1). The results were 
found to superpose very well with those shown m Fig. 3. This check ensured that the 
code was (relatively) free of errors. Fig. 4 shows the 3-dimensional plot characterizing 
the non-dominating Pareto optimal points for Problem 4. This surface has been obtained 
by curve-fitting the several Pareto optimal points actually obtained. It is observed that 
optimal solutions with ppm levels of CO m the flue gas are associated with high air flow 
rates, and yield larger amounts of gasoline. In contrast, optimal solutions with higher (~7 
- 8 %) amounts of CO in the flue gas are associated with lower air flow rates. This is 
similar to what was observed for Problems 1 and 2. Interestingly, it is found that both the 
complete-combustion and the partial-combustion zones are obtained in the present 
problem, simultaneously. It can be observed from Fig. 4 (and detailed tabulated 
information; not being provided here) that every point on the Pareto surface is equally 
good (non-dominating) since any point may be worse than another in terms of one or two 
of the three objectives, but is superior in terms of least one other objective function. If we 
look at the surfe.ce plot closely, we see several maxima and a few adjacent minima. 
Table 6 shows details of two sets of such maxima-minima combinations. Set I in this 
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Figure 4 Suifece plot for 3-D problem (Problem No. 4) 
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TABLE 6 

Details of Two Sets of Chromosomes in Problem 4 



Set 1 

Set 2 

% CO in the flue gas 

0.896 

0.897 

4.268 

5.212 

Gasoline yield (%) 

42.62 

35.62 

41.49 

32.31 

Air flow rate (kg/s) 

20.92 

19.42 

18.56 

17.33 

Air preheat temperature (K) 

456.45 

471.99 

453.01 

458.72 

Feed preheat temperature (K) 

666.75 

596.73 

665.45 

584.84 

Catalyst flow rate (kg/s) 

285.04 

202.93 

261.77 

192.83 




table shows that the air flow rate is lower (better) when the gasoline yield is lower 
(worse), for almost the same amount of CO in the flue gas, i.e., the peak is superior in 
terms of the gasoline yield, but the minima is better in terms of the other objective 
function (air flow rate). A small decrease in the air flow rate (from 20.92 to 19.42 kg/s) is 
associated with a sharp decrease in the catalyst flow rate (from 285.04 to 202.93 kg/s), 
which is responsible for the sharp decrease in the gasoline yield (from 42.62 to 35.62 %). 
Set II in Table 6 shows similar results. This is because of the high sensitivity of the model 
to the catalyst flow rate in certain regions of the parameter-space. Similar sensitivities 
were reported by Arbel et al.^'^ It is difficult to ascertain, at this stage, if this high 
sensitivity is real or just a weakness of the model being used here. An improved 
regenerator model along with the 10- or 12-lump kinetic scheme for the riser will, 
perhaps, lead to an elimination of such sharp negative peaks in the multiobjective 
optimization functions. 

Some regions of the ridge in-between nearby peaks are shown by rectangles in 
Fig. 4. There are no optimal points in these regions, at least from among the 100 
chromosomes obtained in this study. It must also be mentioned here that the computed 
Pareto solutions occur only in the inner zone in Fig. 4 (where all the maxima and minima 
lie), and that no solutions exist in the adjacent flat surfeces in Fig. 4 (these are artifacts 
generated by the graphical package). In other words, in Fig. 4, the Pareto points 
essentially form a series of several hillocks of decreasing heights (with a few nearby 
minima), starting from the low CO-high air flow end (region A) and curving gently 
towards the high CO-low air flow end (region B). 

If we assign some extra importance to the gasoline yield (a decision-maker’s 
post-Pareto job) than used for generating Fig. 4, the several hills in this plot would 



become a set of several preferred solutions. The details corresponding to several of the 
peaks in Fig. 4 (some indicated by +), have been plotted as simpler 2D plots in Fig. 5. 
The two objective functions (gasoline yield and air flow rate) for these selected peaks are 
shown as a function of the third objective function (% CO in the flue gas) in Figs. 5 a and 
b, respectively. All the decision variables corresponding to these peaks are plotted in 
Figs. 5 c-e against percent CO in the flue gas. Several important state variables associated 
with these points are shown in Figs. 5 f-q. It is observed that for the optimal solutions 
shown in Figs. 5 a and b, the catalyst flow rate and the feed preheat temperature lie near 
their upper bounds. This is because high (~ 40 %) values of the gasoline yield require 
high values of both these variables. The constraint on Cigc is taken care of by the air flow 
rate. The constraint on Trgn can be satisfied either by adjusting the air flow rate or by 
adjusting the air preheat temperature. Minimization of the air flow rate leaves Tat as the 
only variable to satisfy the constraint on Trgn. The plot of the air preheat temperature (Fig. 
5 c) shows some scatter. The trends of the total conversion, coke, dry gas and LPG yields 
and the riser bottom temperature (Figs. 5 g, f, h, j and k, respectively) are similar to those 
shown by the gasoline yield. As the percent CO in the flue gas decreases (moving jfrom 
the partial-combustion zone to the complete-combustion zone), the ah flow rate is found 
to increase. This leads to an increase in the oxygen available at the end of the dense bed 
and its subsequent use in the after-burning of CO. This simultaneously results in an 
increase in the temperature rise in the dilute phase of the regenerator (Fig. 5 i), as well as 
a decrease in the coke on the regenerated catalyst (Fig. 5 1). Fig. 5 can be used to 
advantage by a decision-maker to select an appropriate ‘preferred solution (point of 
operation), based on his or her intuition, something that is often not quantifiable. 
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Figure 5 (contd.) Pareto set, decision variables and state variables of selected optimal solutions 
(peaks) of Problem 4 (3D problem) 




Appendix V shows the change in the objective functions over the generations for the 
Problem 4. 

We also studied the effect of varying the computational parameters'*^ associated 
with GA on the 3-D Pareto set (results not shown for the sake of brevity). It is observed 
that by changing the mutation probability, Pm, from 0.05 to 0.025 (refer Appendix VI), 
there is hardly any effect on the Pareto-optimal solutions (the Paretos superpose). 
Similarly, changing the crossover probability. Pc, from 0.95 to 0.975 does not have much 
effect on the scatter of the optimal Pareto-solutions (reger Appendix VII). In the first five 
generations (Problem 4), the scatter in the Pareto set reduces as the generations evolve. 
However, after this relatively small changes occur. 



CHAPTER 4 


CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 

The present study provides a set of optimal solutions for FCCU operation in an existing 
industrial refinery. The procedures are quite general and can be applied to any other unit. 
The optimal results can be broadly classified into two regions, the full-combustion and 
the partial-combustion modes. In order to study both of these simultaneously one has to 
solve a three-objective function optimization problem. This study provides a thorough 
understanding of all aspects of FCCU operating under optimal conditions. Pareto 
solutions provide an opportunity to an engineer to operate an industrial unit optimally. 
Recommendations for Future Work 

The present work should further be extended in order to study the 
multiplicities observed in the 3-D problem in detail. The present off-line optimization 
study should further be extended to on-line optimization. A more rigorous model should 
be used for the prediction of all the unaccounted variables in the present study. 
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(r; as in Eqs. A14 - A16, but with appropriate parameters/variables) 
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Appendix II 


Elitist Non-dominated Sorting Genetic Algorithm (NSGA-II) 

□ Main criticisms of NSGA: 

> Computational complexity of sorting 

> Specifying sharing parameters 

> Lack of Elitism (best parents may get lost) 

• Total complexity of NSGA-II: M(N**2) 

• Total complexity of NSGA: M(N**3) 

M: number of objectives 
N: population size 

1 . A fest non-dominated sorting approach in NS G A-II : 

In order to sort a population of size N depending upon the principle of non- 
domination, each chromosome is compared with aU the other chromosomes to find if 
it is dominated. 

Non-domination criteria: 

A chromosome is better than any other chromosome if at least any one of its fitness 
values is better than the corresponding fitness value of the other chromosome and no 
fitness value of the first chromosome is worse than any of the corresponding fitness 
value of the second chromosome. 

This approach in NSGA-II is similar to that of NSGA. Only difference is a better 
book keeping strategy is used to make the algorithm faster. In this approach, every 
chromosome is compared with partially filled population for non-domination in 


lii 



contrast to what we have in NSGA where it is compared with all the other 
chromosomes. We temporarily put one chromosome (p) in a separate box, and then 
we add one chromosome (q) from the population to that box, and compare q with all 
members in the box for non-domination. If any member in the box is dominating q, q 
is removed. If q dominates any member in the box, that member is removed. But if q 
is non-dominating to all members in the box, then q is put in the same box and 
continues -with the other members. This is how the temporary box grows with non- 
dominated chromosomes. When all chromosomes are checked, the chromosomes left 
in box are assigned with rank, Irank (Here Front 1 is created so rank assigned is also 
1). All the other members (not present in Front 1) are used for evaluating other non- 
dominated fronts and assigned subsequent ranlcs. 

2. Density Estimation: 

One of the criteria for judging the betterment of any evolutionary algorithm is how 
much diversity it can preserve among its solutions. The original NSGA used the 
famous sharing function approach for preserving diversity among its population. 
Major difficulties with sharing approach are its performance depends upon the ctshare 
value and the overall complexity of this approach is quite higher. In NSGA-II, the 
sharing fiinction approach is replaced with a crowded comparison approach in order 
to overcome all the above-mentioned difficulties. 

The density of each chromosome is estimated by calculating crowding 
distance. In this, all the chromosomes belonging to a particular front (or all those 
having the same rank) are sorted according to each objective function value in their 
ascending order of magnitude. Boundary solutions are assigned with am infinite 



distance value (practically speaking a very large value). The intermediate solutions 
are assigned with a distance value equal to the absolute difference in the function 
values of two adjacent solutions. This is done with all the other objective functions. 
The overall crowding distance (Idist) is calculated as the sum of individual distance 
values corresponding to each objective. This method is continued for all the other 
fronts. 

3. Initially (at generation = 0) all the parent solutions are sorted based on the non- 
domination. Each solution is given a rank based on its non-domination level. All 
these solutions undergo binary tournament selection (Binary tournament selection is 
carried out using crowded comparison operator, see stage 7), crossover and mutation 
operations and produce child population, which is equal in size to that of parent size. 

4. In NSGA-II, parent solutions and child solutions are mixed together to ensure elitism. 
Elitism is a concept that helps in preserving the elite members of the population 
throughout all the generations. By mixing parents and children, elitism is ensured. 

5. Once the elitism is done, the whole population (now of size 2N) is sorted on the basis 
of non-domination approach as explained in stage 1. Solutions belonging to the first 
non-dominated front (or rank) are the best solutions and should be preserved. 
Crowding distance is calculated for all the solutions belonging to different ranks as 
explained in stage 2. 

6. Now, if the number of solutions belonging to first non-dominated front (rank 1) is 
smaller than the total population size (N), all the solutions are included in the new 
population. The remaining members for the new population are selected from 
subsequent non-dominated fronts in order of their ranking. This is continued till we 



fill up the new population with the members equal to the original population size (N). 
To ensure exactly N members in the new population, the last front is sorted on the 
basis of crowded comparison operator as explained in stage 7. 

7. Crowded comparison operator: 

This is used to compare the two solutions at any stage in the algorithm. We know that 
every solution has non-domination rank and crowding distance associated with it. The 
operator is defined as; 

If we have two solutions with different ranks, we select the solution with the lower 
rank (means the one from the better non-domination front) 

If we have two solutions belonging to same non-domination front, we select the one, 
which is located in a lower crowded region. 

It can be seen from the algorithm that the sorting is done with a better and 
faster approach as compared to original NSGA. The diversity is preserved using the 
crowded comparison operator instead of using the sharing approach as used in 
NSGA. There is no need to specify any sharing parameter. Even though the crowding 
distance is calculated in the objective function space, it can also be calculated using 
the parameter space. 



Flowchart: 












Appendix III 


Comparison of Model Results with Industrial Data 



Data Set 1 

Data Set 2 

Data obtained 
from Industry 

Results obtained 
After tuning 

Data obtained 
from Industry 

Results obtained 
After tuning 

Feed flow rate 
(kg/s) 

"28.33 

28.33 

29.08 

29.08 

Catalyst flow rate 
(kg/s) 

208.33 

208.33 

205.00 

205.00 

Feed preheat 
temp (K) 

617.4 

617.4 

620.7 

620.7 

Riser pressure 
(atm) 

2.457 

2.457 

2.506 

2.506 

Air flow rate 
(kmol/s) 

0.559 

0.559 

0.5711 

0.5711 

Air preheat 
temp(K) 

490.3 

490.3 

493.3 

493.3 

Regenerator 
pressure (atm) 

2.588 

2.588 

2.638 

2.638 

Riser top 
temperature (K) 

765.5 

763.0 

766.6 

767.27 

Regenerator 

temp(K) 

930.2 

925.0 

934.6 

935.79 

Flue gastemp(K) 

945.4 

945.34 

946.8 

953.38 

Flue gas 
oxygen (%) 

0.00029 

1.21e-l 

0.0002 

5e-2 

Dry gas % 

3.1 

3.6 

3.8 

3.635 

LPG % 

12.1 

12.66 

12.6 

12.821 

Gasoline% 

32.6 

34.92 

35.2 

35.454 

Gas oil % 

48.1 

44.287 

44.1 

43.594 

Coke % 

4.1 

4.519 

4.3 

4.4932 







Data Set 3 

1 Data Set 4 

Data obtained 
from Industry 

Results obtained 
After tuning 

Data obtained 
from Industry 

Results obtained 
After tuning 

Feed flow rate 
(kg/s) 

29.08 

29.08 

29.11 

29.11 

Catalyst flow rate 
(kg/s) 

208.33 

208.33 

208.33 

208.33 

Feed preheat 
temp (K) 

621.4 

621.4 

625.1 

625.1 

Riser pressure 
(atm) 

2.545 

2.545 

2.546 

2.546 

Air flow rate 
(kmol/s) 

0.5759 

0.5759 

0.5707 

0.5707 

Air preheat 
temp(K) 

494.4 

494.4 

493.9 

493.9 

Regenerator 
pressure (atm) 

2.678 

2.678 

2.68 

2.68 

Riser top 
temperature (K) 

768.5 

769.83 

768.8 

770.139 

Regenerator 

temp(K) 

937.5 

937.37 

937.5 

936.88 

Flue gastemp(K) 

949.1 

952.0 

948.2 

948.4 

Flue gas 
oxygen(%) 

0.00028 

3.92e-2 

0.00031 

3.45e-2 

immmm 

3.8 

3.7 

3.4 

3.69 

LPG % 

12.6 

13.0 

12.4 

13.05 

Gasoline% 

35.2 

36.1477 

34.0 

36.13 

Gas oil % 

44.1 

42.545 

45.6 

42.594 

Coke % 

4.3 

4.543 

4.6 

4.522 




conversion (%) feed preheat temperature (K) 


Appeudix IV 


Pareto set, decision variables and state variables for optimal solutions 
of Problem 3 

population size 50; number of generations 50; mutation probability 0.05 
crossover probability 0.95 
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Appendix V 


Evolution of Pareto (Problem 4) over generations 
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Appendix VII 


Effect of crossover probability on Pareto set, decison variables and state variables 
of Problem 4 (population size 50; number of generations 50; mutation probability 0 


P = 0.95 (filled circles); P = 0.975 (unfilled circles) 
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